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Abstract- The digital frequency synthesizer (DFS) is an integral part of the 
programmable local oscillator (PLO) which is being developed for the 
NASA's Deep Space Network (DSN) and radar astronomy. In this report, 
the theory of operation and the design of the DFS are discussed, and the 
design parameters in application for the Goldstone Solar System Radar 
(GSSR) are specified. The spectral purity of the DFS is evaluated by 
analytically evaluating the output spectrum of the DFS. A novel 
architecture is proposed for the design of the DFS with a frequency 
resolution of 1/2^8 of the clock frequency (0.35 fiHz at 100 MHz), a phase 
resolution of 0.0056 degrees (16 bits), and a frequency spur attenuation of 


-96 dBc. 
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L Introduction 


This report describes the theory and the design of a digital frequency 
synthesizer (DFS) used in the programmable local oscillator (PLO) for 
applications in the NASA’s Deep Space Network (DSN) and the Goldstone 
Solar System Radar (GSSR) [2,3,4]. The DFS provides a low cost solution for 
frequency synthesis with the following features: 

•Controllable frequency and phase with extremely high resolution 
•Long term stability in the output frequency 
•Low phase noise and high spectral purity 
•Phase continuous frequency adjustment 

In Figure (1), the overall configuration of the X-band exciter proposed in [1] 
is shown. This exciter is used for transmission of a coherent X-band signal 
from the DSN station, to a distant solar system planet. The reflected signal 
from the planet is used for generating radar images of the planet [2,3], 


O 

%% 
iS *-> 

d? S3 
Pm m 

% PP 

e E 

94 W 


S’ 



Figure (1). X-Band Exciter for the Goldstone Solar System Radar 

Transmitter 

The design of the exciter is based around using a high-resolution PLO with 
controllable phase and frequency. The output of the PLO is ideally a single 
carrier with a frequency range of 10-20 MHz. This signal is used as a 
reference by an HP-86621 synthesizer; the output of the synthesizer is up- 
converted with an ultra-stable m^ser clock provided by the station. This 
signal is then modulated and transmitted in the X-band. 
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The main focus of this work is the design and analysis of the PLO shown in 
Figure (1). In Figure (2), the block diagram of the PLO is depicted. The 
PLO is composed of a DFS and digital-to-analog conversion module, and is 
controlled by a host via an HP I 11 Interface Bus (HPIB). 
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Figure (2). Programmable Local Oscillator Block Diagram 

The PLO controller software controls both the frequency and the phase 
from the PLO, via the HPIB. Other control functions of the host via the 
HPIB are described in more detail in [4]. 

II Organization of This Report 

In sections II through IV, the design and analysis of the DFS are covered. 
In section V, the design requirement in application for GSSR is outlined. 
In section VI, it is shown that none of the off-the-shelf DFS modules 
available in the commercial market satisfy our design requirement. In 
section VII, an architecture for the DFS is proposed that satisfies our 
requirement. 


II Digital Frequency S ynthesizer 

With the advent of high-speed digital signal processing hardware, today 
frequency synthesis for intermediate frequencies (IF) can be achieved in the 
digital domain. A DFS uses a single reference frequency to generate a 
range of output frequencies. In general, an ideal frequency synthesizer 
would generate a single sinusoid at a desired output frequency denoted as 
S(t) =Sin(2nF 0 t) . The underlying theoretical model of the DFS is shown in 
Figure (3). 


I HP™ is the trademark of the Hewlett-Packard Co. 
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Figure (3). Theoretical Model for the Digital Frequency Synthesizer 


The phase and frequency of the output sine wave are controlled by the 
frequency control word, denoted as F r . This value is used by the phase 
generator to compute the phase sequence (Kn ); this sequence is then 
mapped around the quantized unit circle with each clock pulse, which 
occurs at frequency F c . The granularity of the phase space depends upon 
the number of quantization intervals around the unit circle, and is taken to 
be equal to 2^ in Figure (3). Ideally, the sequence S(n) is the sampled 
values of S(t) at the discrete sampling points nT, where each sample is 
taken at the rate of T=1 /F c , i.e., S(n)=Sin(2riF 0 nT). In the next section, the 
relation among F c> F 0 and F r is clarified. 

m. Theory of Operation 

The implementation of the DFS shown in Figure (4) was first introduced in 
[8]. This approach employs a binary accumulator for the phase generator, 
and a look-up table for converting 6(n) to S(n) . 



Figure (4). Implementation Model for DFS 


The theory of operation for this model is as follows: an L-bit binary adder is 
incremented using modulo 2 L addition by the frequency control word: F T , 
and the accumulated result is stored in the phase register as %(n). The 
sequence £(n) is truncated from L-bit to W-bit to produce 6(n), which is used 
to address a look-up table stored in a read-only memory (ROM) to output the 
desired output samples S(n). In order to provide sufficient resolution in 
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most applications, the binary accumulator size must be in excess of 20 bits 
long. It is practically impossible to use a ROM with a depth of 2 20 for 
converting the phase sequence £(n) to the sine samples. Hence, it becomes 
necessary to truncate the output of the phase generator by B-bit to a 
practical range W=L-B. The truncation £(n) to (Kn), is equivalent to a 
modulo-2^ s ) division (or shift), as shown in Figure (4). The truncated 
phase sequence (Kn ) is used as an address to the 2 w xD-bit-wide ROM which 
maps (Kn) into the desired sequence S(n). 

The output frequency of the DFS, denoted as F 0 is related to F c , and F r by 
considering the recursion relation for the phase sequence <fj(rc), which is 
£(/i+l)=(£(n)+F r ) modulo-2^. The rate of change in each unit time for <f;(n) is 


A Un) 

At 


= FF 


( 1 ) 


In case of infinite precision in the ROM and no phase quantization with 
( Kn)=^(n ), the output is S(n) -Sin(2nF 0 nT) , where the output frequency F 0 
for the sine function is 


0) 

o 

2k 


MM 
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( 2 ) 


Note that F x is an integer valued variable, and F 0 is a real valued variable, 
thus for a given fixed clock frequency and accumulator length L, the 
frequency control word is FV=L2AF 0 !F c J, where tad denotes the floor of x. 

From the sampling theorem, it is well known that the maximum frequency 
of a sampled sinusoid is bounded by F c f 2. Hence, we have 

F 0 Uax = ^ => from (1) 1 < F r < 2 t_1 (3) 

From equation (3), it is concluded that it is possible to generate 2^ _1 ) 
possible output frequencies with a frequency resolution of F C /2 L . 

To explore the spectrum of this sequence consider the sequence £(n) which 
is periodic. The period here is defined as the smallest integer N such that 
%(n)=%(n+N). This period is N=2 L /gcd(F T £ L ), where gcd(a,P) denotes the 
greatest common divisor of integers a and p. Thus, N will be equal to the 
period of the output sine wave if and only if F r =gcd(F r ^) , or equivalently 
F r is purely a power of two. Since %(n) is periodic with period N, it can be 
represented by a Fourier series composed of N-harmonics with the first 
fundamental frequency at F r /2^. We refer to these N - 1 harmonics 
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(excluding the first fundamental frequency) as "frequency spurs", which 
form the induced interference with the desired output signal. 

In case W<L, the truncated phase sequence 0(n)=[£(n)/ where [jt/y] 
denotes taking the integer part of the division of x by y, has the same period 
N as £,{n). Therefore, the sequence (Kn ) can also be represented by a 
discrete spectrum of a periodic sequence consisting of N spectral lines. 
Furthermore, if F r and 2 L are mutually prime (i.e., F T is odd), the number 
of spectral lines is exactly 2 L , which in practice is a very large number. 

The size of the arithmetic field used to compute the phase sequence (Kn ) 
directly impacts the frequency spur characteristics of the DFS. When using 
the binary field to compute this sequence, the DFS's output sine wave 
consists of N=2 L /gcd(F r ,2 L ) harmonics with the first fundamental 
frequency at FfF^L, and the number of possible output frequencies is also 
always restricted to be powers of two, if binary number representation is 
used for representing (Kn). It should be noted that this restriction can be 
avoided by using other arithmetic fields such as the decimal field, or 
alternatively other forms of number representation such as binary coded 
decimals (BCD). 

m.l Output of DFS D ue to Phase Truncation 

For simplicity assume without loss of generality that F c = 1. When the phase 
sequence Efn) is truncated to W-bit, let L=W+B, then the the output 
sequence is 


S(n) = 5m 


2k- 


Z(n) 


T\ 




(4) 


The integer part of the ratio in the numerator of (4) can be written as 


2 b 



+ £(n) 


(5) 


where e(/i) is an error sequence corresponding to the remainder of the 
expression in (5), and its value is always bounded by one. Using equation (5) 
in (4), the output sequence can be expanded as 
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Sin 
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400 ' 
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= Sin (j£T Z in) Cosine (n ) j + Cos (i^r £ ( n )j Sin {^r £ ( n )j (6) 


Under reasonable circumstances, 2 L »1, the magnitude of the sequence 
e(n) is small so that Sin(2nein)/2^) will be small compared to 
Sin(2n^(n)/2 L ) and Cos(2nein)/2 w ) is close to 1. Hence, we can make a first 
order approximation to the output sequence S(n), denoted as S\(n) to be 


S,(/i) = Sin 


r 2n%(n) 

2 L 


\ 

) 


(7) 


which is the desired output sequence. To estimate the frequency and the 
size of the largest spur which is denoted by y/(n), we look at the next 
approximation, and note that 



Since e(n) is small, Sin 2 (nein)/2 w )«Sin(2n£(n)/2^) for n such 
Sin(2n^(n)/2^ , )<Cos(2n^(n)/2 L ), the size of the largest spur, by 
approximating the second term in (8), and bounding the magnitude of this 
term as follows 


Sin^£(n)^J 


< 2££(rt) < 


n 




.W-l 


(9) 


The first inequality is valid since 2~ w e(n) is small, and the second 
inequality is due to the fact that e(n ) is always bounded by one. These 
approximations are not valid for certain frequency ranges, and therefore it 
becomes necessary to analytically evaluate the output spectrum of DFS. 

m.2 Output Spectrum of DFS 

To analyze the spectral purity of the DFS, it is desirable to derive an 
analytical expression for the size of the spurs as a function of frequency. In 
reference [5], the authors attempt to drive an analytical expression for the 
output of DFS, however, their approach has a number of flaws [6]. 


DFS for Radar Astronomy 


Jet Propulsion Laboratory 


It is shown in Appendix A that the ratio of the power of the output frequency 
to the power of the largest spur denoted as p is 

p dB = 20 Log , o (2 W ) (10) 

This equation is shown to be valid only for odd values of F r , and the location 
of the largest spurs are m such that m /F r =+1+2 W , and m /F r =-1+2 W . As it 
will become evident later, the requirement of F r being odd is not restrictive, 
since a small modification in the architecture of DFS could assure this 
condition without any impact on the performance, resolution or the spectral 
purity. 

IV. Simulation Results 

In cases when L is a small number (<16) it is possible to obtain the output 
spectrum of DFS by simulation. The DFS shown in Figure (4), was 
simulated on Mathematica® 1 [9], which is equipped with Fourier 
transform routines. The output of this program is shown in Figure (5) for 
the parameters are L=12, F c = 1 and W= 5 . This example was chosen to 
verify the result against those of reference [5], which was obtained directly 
from the output hardware of a DFS. 



1 Mathematica® is a trademark of Wolfram Research Inc. 
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IV. 1 Effect s of Windowing 

The output sequence from the D FS is periodic with period N. When 
simulating the DFS, it is not practical to generate and store the whole 
sequence of length N, since the range of values for N in our case of study in 
some cases is in excess of 2 32 . Furthermore it is not feasible to compute the 
discrete Fourier series for very large N. Thus, it becomes essential to 
window the output sequence S(n) to a practical range. 

Let S~(n ) denote the windowed sequence, .i.e., 

S(n) = ^S(i)w(n-i) (11) 


Consider a window of length M , that is 

r -m ' m 

w(n) - 1 * ’ 2 ~ H ~ 2 (12) 

[ 0, Otherwise. 

Let <&i(o)=j(S(n)) and <P~(Q))=j(S~(n)) each respectively represent the Fourier 
transforms of S(n) and S~(n), where y(.) denotes the discrete Fourier 
transform operator. Since <Kco) is periodic with period N, it can be written 
as 

yV-1 

(P(0)) = (13) 

, -N - 1 N 


Here S(.) denotes the delta-Dirac function, and the Fourier series coefficient 
ak is 


N- 1 

2 .2nnk 

%S(n)e~ (14) 

N-\ 

2 


The Fourier transform of the windowed sequence can be evaluated as 
<jy-(G))= &(cd)®W(G)), where ® denotes the convolution operation, which is 

K 

O(o>)= j^W(o)-y)a k 5(y-co k )dy (15) 

-* k 

By changing the order of summation and integration in (15) and evaluating 
the output spectrum at each point a>i=27d/M > with i=-MI 2, (M+l)/2 ,..., 
0,...,(Af-l)/2,M/2 (15) becomes 


DFS for Radar Astronomy 


Jet Propulsion Laboratory 


III'! IT I 1 IP! I 



O(oo)= 2, a k W(a)- : 

, 1W-D 


<$^(< 0 ) in (16) is a distorted version of <Hco) due to the spectral leakage caused 
by windowing. To illustrate the effect of spectral leakage consider the 
Fourier transform of a rectangular window of length M which is WXco) 
=Sin(oMI2)l Sinioj! 2) [10]. The amplitude spectrum ofW(co) is shown in 
Figure (6. a). The zero crossing of this function is at 2jd/M. It is evident 
from equation (16) that the output spectrum is the superposition of the 
original spectrum and the basis function Sin{(£MI2)ISin{(ol2) as shown in 
Figure (6.a) shifted at each point 2 ni/M, with i=-M/ 2,-(M+l)/2,...,0,...,(Af- 
l)/2,M/2 , as shown in Figure (6.b). 


I Bran 










l ?=sk i^viiiiS i 


Figure (6). Fourier Transform of the a* Window b. Basis Function c. 

In Figure (6.c), the sampled spectrum is shown for every 4-th impulse. At 
each discrete point spaced at 2n!M , the original spectrum is superimposed 
with all the other (M-l) sidelobes. In Figure (7.a and b), the distorted 
version of the output spectrum is shown for the example of Figure (5), 
where AT=2048, and M= 256. 
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Figure (7). Effect of Windowing on the Output Spectrum for cases B. and C. 

of Figure (5). 


The spectral leakage causes a bias in the amplitude and the location of the 
harmonics. To reduce this effect, the window jshould exhibit low-amplitude 
sidelobes far from the central lobe, and the transition of to the low sidelobes 
should be very rapid. Hence, a design criterion to suppress the leakage is 
the peak sidelobe level, relative to the main lobe. For the rectangular 
window shown in Figure (6. a), the highest sidelobe level for the rectangular 
window is -13 dB. To assess the spectral purity of the DFS in our 
application, it is necessary for the highest sidelobe level to be around -96 
dB. A class of windows known as the Kaiser-Bessel window exhibit this 
desired level of sidelobe attenuation. This window is defined by 




w(n) = 



/<,(*«) 


0 <«< 


M 


where, /„(*) = J 


2 * 


£S(2*9 


2 ' 


(17) 


The parameter na is half of the time bandwidth product, and with a=4.0 
the highest side-lobe level is -97 dB. In Figure (8) both the window sequence 
and its amplitude spectrum are shown when a=4.0. 



Figure (8). Kaiser-Bessel Window ct=4.0 
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V. PFS Requirement for Our Case of Study 

The design requirement for the DFS is summarized in Table [1]. 


Desicm Reauii 

rements 

Output Frequency 
Frequency Resolution 
Settling Time Maximum 
Spectral Purity 
Dynamic Range 
Spectral Continuity 
Phase Resolution 
Phase Stability 

10-20 MHz 
1 pHz 
5 msec 
-70 dB/c 
-72 dB 

Phase continuous 
0.2° 

l°/24 hours 


Table [1]. Design Requirement for GSSR 


This requirement is extracted from references [1,4], where the advanced 
exciter is described and the advantages of this requirement for the DSN and 
radar astronomy are discussed. 


VI. Evaluation of Off-the Shelf DFS Modules 

At the time when this study was undertaken the two fastest DFS modules 
available off-the-shelf were Qualcomm 1 Q2334- and Sciteq 2 VDS-30000. 
However, it works out that neither one of these products satisfies our design 
requirements as it is described below. 

Qualcomm Q2334: 

This single chip also has two problems: 

•It only provides a frequency resolution of 1/2^2x30x1 0^=7mHz, which is not 
sufficient for our application. 

•The maximum clock frequency is 30 MHz, which is not sufficient for our 
application. 

Sciteq VDS-3000; 

This module provides the required frequency resolution with (from 
equation(2)) l/2 48 x64xl0 6 =0.2pHz which is sufficient for our application. 
However, there are two problems with this module: 

•The phase resolution is only 1.4° (360/2 8 ), which is not sufficient for our 
application. 


1 Qualcomm® is a registered trademark for Qalcomm Corp. 

2 Scitec® is a registered trademark of Scitec Corp. 
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•The dynamic range of the output is only 8 bits, which is also not sufficient 
for our application. 

However, the phase resolution which is 360/2 16 = 0.00549316° and the 
dynamic range, which is 12 bits, are appropriate for our application. 


VII.DFS Im plementation 

In this section, an architecture is given that satisfies our design 
requirement, and it is also well suited for VLSI implementation. This 
architecture employs standard off-the-shelf Bipolar ECL (Emitter Coupled 
Logic) components for implementing the DFS model shown in Figure (4). 
The proposed architecture has the features and specifications shown in 
Table [2J: ; 


f Design Features 

Output Frequency 

0-50 MHz 

Frequency Resolution 

0.35 (iHz at 100 MHz Sampling Rate 

Spectral Purity 

-96 dB/c Worst Case 

Dynamic Range 

Min -72dB Max -96dB 

Spectral Continuity 

Phase continuous 

Phase Resolution 

0.0056° 

Hardware P 

requirement 

Accumulator Length 

48 bits 

Sin (.) Look-up Table Resolution 

16 bits 

Maximum Clock Frequency 

110 MHz 

Output Resolution 

12(16 bit 


Table [2]. Design Features and Hardware Requirements 

The resulting architecture not only satisfies our design requirement, it is 
also the fastest DFS with highest resolution ever built, thus providing an 
effective approach well into the next century. 

VII I Spectral Purity 

In this section the Kaiser-Bessel window outlined in Section IV. 1 is 
employed to simulate the output spectrum of the DFS. In Figure (9), the 
simulated spectrums for our design parameters (F c = 50 MHz, L=48, 1^=16, 
M= 1024) for two different frequencies are shown. In Figure (9.b), the output 
frequency was chosen so the period of the output sequence is 2^. 
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c 

Figure (9). Output Spectrum for Our Case of Study 


As before, we conclude that the spectral purity of the DFS is dependent on 
the frequency control word F r . From this simulation and equation (12), it is 
concluded that the spurs are well below our design requirement of -96 dBc. 

Vn.2 DFS Architecture 


Our proposed architecture for the DFS is shown in Figure (10). The design 
methodology is a fully synchronous sequentially pipelined architecture, 
synchronized with the system clock. The main difference between this 
architecture and the original model shown in Figure (4) lies in the manner 
in which the Sin(.) component of the phase accumulator is computed. The 
phase generator is a standard 48 bit binary accumulator. Frequency control 
is provided by loading a 48 bit value computed from equation (2). Phase 
control is achieved by preloading the upper 16 bit of the 48 bit internal 
register of the accumulator. Once the device is enabled the content of the 
internal register is continuously incremented by the value of the frequency 
control register. The most significant word (16 bits) of the internal register 
denoted by the index i is used to address the look-up table which output is 
Sin(i2ji/2 16 ). 
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Figure (10). DFS Architecture 

One of the fundamental difficulties in designing the DFS is the depth of the 
look-up table to generate the sine samples of the phase sequence. A simple 
method to reduce the table size is to decompose the output of the phase 
generator in binary representation, where 

9 = 2 15 d 15 + 2 14 d 14 -t-’-+2 1 d 1 + d Q where, d i e {0,1} V,. (18) 
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Define a and P in terms of 6 as 


a = 2 l5 d x5 + 2 14 d, 4 + 2 n d n + 2 n d X2 + 2 11 d n + 2 i0 d x0 + 2 9 d 9 + 2 % d % 
p = 2 1 d 1 +2 6 d e +2 5 d 5 + 2 *d A + 2 V 3 + 2 2 d 2 + 2d x +d 0 (19) 


Using this decomposition we have 


Sin(6) = Sin( 
Sin( 


a+p 

2 16 


/" ->7 


2tt) = Sin 


2 7 djj +• • • +d% 2 1 cLj+‘-'+d 0 


2^ 


,16 


)2n 


( 20 ) 


• /2 7 d, 5 +” - +do ^ / 2 7 d 7 H — i-d 0 „ . _ w 15 

’ " ; 2 n)Cos ( — - z —^ 2?r) + Cost 


2’d.+-+j, 2 k ) 


The computation of Sin(i2n/2 ie ) is now reduced to two look-up tables with a 
depth of 256 points (as opposed to 65536 points), with an additional 
complexity of two multiplications and an addition. Equation (20) is 
essentially equivalent to quantizing the unit circle into two levels. The first 
level is coarse quantization with a quantization interval of 2rc/256 as shown 
by 'x'-mark in Figure (11), and the second level is fine quantization with a 
quantization interval of 2 tc/65536. Li Figure (11), the fine quantization 
intervals are shown by solid lines. In the block diagram of Figure (10), the 
in-phase and quadrature components of the fine quantization points are 
denoted as SinF and CosF correspondingly, and simply as Sin and Cos for 
for the coarse quantization points. 
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To compute equation (20), it is necessary to perform two multiplication and 
a single addition, organized as a sequential pipeline. This is done by 
multiplying one component from each set of fine and coarse angles, as 
shown in Figure (11), resulting into a 32-bit value in which the higher 16- 
bit values for in-phase and quadrature (I&Q) components are added to yield 
the desired result. The multiplication operation is performed using a 
standard off-the-shelf multiplier, with the rounding feature enabled at the 
15 th -bit. The rounding from the least significant bit is also performed at the 
last stage as a part of the full adder. 

To assure the spectral purity of the DFS expressed in terms of p in equation 
(12), the following minor modification to the control word register can be 
made to force the accumulator to always make odd frequency steps. 


LSB from the 
Control Word 



Figure (12). Modification to the Content of Phase Generator to Force F r to 

Become Odd 

This modification amounts to using the least significant bit of the phase 
accumulator, as if the word length of the accumulator is L+l and the least 
significant bit is always set to one. This logic circuit essentially doubles the 
frequency resolution of the phase accumulator by adding an extra bit to its 
word length, and forcing the input frequency control word F r to be used 
internally as if it is an odd number. This is performed without any 
compromise to the frequency resolution of the DFS. The input to the carry- 
in of the 48-bit accumulator is toggled from one to zero when the least 
significant bit of the phase accumulator is one, and is toggled from zero to 
one when the least significant bit of the phase accumulator is zero. This 
operation is effectively equivalent to alternately flipping the carry-in to the 
accumulator in each clock cycle. It is pointed out in reference [5] that this 
circuitry improves the spectral purity of the DFS 4 dB; however, their 
suggestions are based on some flawed results [6]. Our motivation to use 
this modification is based on our analysis on the spectral purity discussed 
in section III. 

This completes our discussion of the overall architecture of the DFS. In 
the next section, we outline the detail design of the look-up table to generate 
the in-phase and quadrature components. 

Vn.2.a Design of the L ook-Up Table 

The look-up table is composed of the I&Q components of the fine and coarse 
quantization intervals. The coarse quantization intervals span the whole 
unit circle, hence it is only necessary to store a single quadrant for the 
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coarse I&Q values and adjust the sign prior to computing (20). The I&Q 
values for the fine quantization angles are small and, for 16-bit fixed point 
presentation of these values, the in-phase component is between 1 and 
0.99969, thereby requiring only 4 bits of resolution. The quadrature 
component is bounded between 0.0 and 0.0246, thus requiring only 10 bits of 
resolution. The architecture shown in Figure (13) for the look-up table 
takes advantage of these observations. 



Figure (13). In-Phase & Quadrature Component Look-up Table 


In Figure (13), the 16-bit output of the phase generator is split into three 
parts; the least significant byte is used to output the I&Q components of the 
fine quantization angles denoted by a single bus name FSC[.] of a width of 
14 bits which is in a packed format with 4-bit in-phase and 10-bit quadrature 
component. The upper 6-bit is used to address the I&Q components of the 
coarse angles(Sin[], and Cos[]) with 16-bit resolution corresponding to each 
component. 

The format of each output from the look-up table, which forms the input to 
the multiplier accumulator for computing (20), is shown in Figure (14). 
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CosF(.) 


Cos[1 5:0] 



Figure (14). Coarse and Fine I&Q Format 


The format of all the numbers in the look-up table is in two's complement 
format. However for the coarse angles, a scaling factor of (2^ 1 -1)/211 must 
be used to scale all the I&Q values. This is due to the fact the D/A module 
uses offset binary format, and given that this scaling is performed in the 
look-up table, then converting the output sequence from the two's 
complement format to offset binary format amounts to inverting the sign of 
the most significant bit of the results. 


vn.2.b DFS Address Map 

For control operation of the DFS, the following registers are mapped into 
the address space of the controlling host. In Table [3], it is assumed that a 
fixed base address for the DFS is used, which is decoded from the host bus, 
i.e., the addresses are the offsets from the base address and are shown in 
the left-hand column. The second column in Table [3] describes the 
function of accessing each address space, and the last column indicates 
whether the address space is readable (R) or writable (W). 


Add 

Function 

Purpose 

R/W 

0x1 

Fo -Reg 

48 bit Frequency Control Register 

W 

0x2 

Phase Reg 

16-bit Phase Accumulator 

R/W 

0x3 

Out-Reg 

Output Register 16-Bit 

R i 

0x4 

Reset 

Software Reset to zeroize all Regs 

W 

0x5 

Freeze 

Disable Clock 

W 

0x6 

Go 

Enable Clock 

W 

0x6 

Enable 

Device Enable 

W 


Table [3]. DFS Address Map 
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The first three addresses essentially provide the means to control the 
frequency and the phase of the DFS. The last four addresses initiate a 
control function that can be performed under software control from the 
host. The purpose of a software reset is for initialization of the DFS, and 
the freeze and go functions are provided so the DFS can be stopped and go 
anytime during its operation, while maintaining phase and frequency 
continuity. This feature is useful in application of radar astronomy or 
spread spectrum communication. The device enable is provided for 
software control of the device from the host. 


yill. Conclusion 

In the first part of this report, a DFS module which forms the main 
building block for the PLO was analyzed and designed. The spectral purity 
of the DFS was studied using both analytical and simulation methods. An 
architecture for the DFS was outlined which is also well suited for VLSI 
fabrication. This architecture provides a frequency resolution of 1/2 48 of the 
clock frequency (0.35 pHz at 100 MHz), a phase resolution of 0.0056 degrees 
(16 bits), and a frequency spur attenuation of -96 dBc (16-bit output 
frequency). 
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A ppendix A. Deri v ation of Output Spectrum of ] 

Consider the function 


S(n) = Sin 


2n 


(nF) 




where (nF) ^ means we take the representative 




(A.1) 

modulo 2^ of nF among the 


numbers (0,1, ...,2 B -1}. 


Case 1: 

Let F= 1, and consider the function (pin)=[n2^l 2®]. The values of (pin) are 

n< 2 B <p(n )= 0 

2 B <n<2.2 B y(n)=l 


k.2 B <n<(k+l)2 B (p(n)=k . 

So the function <p(n) takes values 0,l,2,...,2W-l. The Fourier transform of 
S(n) is 

S(m) = Y j e 2 * innl2L S(n) (A2) 


Substituting from (A.l) for S(n) we obtain 


2 W -1 ( (*+1)2® -1 \ 
k = 0 V n=k2* 


C- ( lnk \ 

Sln (-^r) 


Using the identity 


(*+l)2*-l \ 

i^ e 2xjm*l2 L 

n=*2* J 


_ gllljmklT 


w 1 — e 




l-e 


2fijm/2 L 


we obtain 


S(m) = 


l-e 

1-e 


2xjm/2 w 2 W -1 


2 nk s 


2n}m\2 L 


. La . 

e SmC-jjr) 


k = 0 


(A3) 


(AA) 


(A. 5) 
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Now 


vT 1 ,2 nk. f27cmk'\ c . (Ink 

e SwC-^r) = 2- Co \”2»” J 5l \^ r 

, Y'V (2nmk'\ (2nk\ 


(A.5) 


Using the two trigonometric identities 


2 CosA SinB= Sin(B+A)-Sin(B-A), 
2 SinA SinB=Cos(A-B)-Cos(A+B), 


we rewrite (A.5) in the form 




It follows that if m*± 1 Mod 2^, then 


S(m) = 0 


(A.7) 


If m =1 Mod 2 W , then 


1 _x,2*/72" 


(A8) 


=-l Mod 2^ then 


S{m) = -2 w - i j 


w - i ,• 1 — e 


-2xjml2 l 


(A. 9) 


Q.E.D 
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Case II; 


F and 2 are relatively prime. The numbers (nF) 2 L are the permutation of the 

numbers {0,1,2 2 L -1} induced by multiplication by F Modulo 2 L . Let us 

write Spin) for S(n) to emphasize the dependence on F. Then Sf(n)=Si(Fn). 
Consequently, 


S F (m) = £e 2 *- /2i S,(n) = W 

n = 0 »=0 


= S l (mF- 1 ) (A 10) 


This settles the case g.c.d{2,F)=l. 


Let us obtain approximations to the largest and the second largest values 
for ISF A (m)l . Clearly the largest value of IS; A (m) I occurs for m=± 1. For 
such values of m we have the approximation 


l-e- 2 * j/2W 

\-e 2K]llL 


2k 



2 l 


(All) 


so that 


I5j(±l)l= 2 L ~ l 


(A 12) 


The second largest value of I S/ A (m) I occurs for m=±l+2 w . Here we obtain 
approximations 


2k 

1 _ g-2*v72 W 

1 _ g 2* y (± l + 2*')/2‘ ~ 2^(±1 + 2' V ) 

2 l 


— * 2 b w therefore 

±1 + 2 W 


S,(±1 + 2 W )| = 2* _i (A 13) 


When F is relatively prime to 2, i.e., g.c.d(F,2)=l, from (A.9), we have 
SF A (fn)=Si A '(mF- 1 ) t hence the largest value of ISjAm) I occurs for 
m=±F Modulo 2 L 


\S f (±F)\~ 2 l ~ j (A 14) 

Similarly, the second largest value of IS A /Km)l occurs for m=F(±l+2 w ), 
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I S F (F(±l + 2 W ))l * 2 b_1 (A15) 

Thus, the ratio p of the power of the output frequency to the power of the 
largest spur which is obtained by dividing (A. 14) by (A. 15), and when 
expressed in dB is 

p=20 Log io (2 W ) (A. 16) 
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